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Abstract 



We address the issue of fluctuations, about an exponential lineshape, in a 
pair of one-dimensional kicked quantum systems exhibiting dynamical local- 
'^ \ ization. An exact renormalization scheme establishes the fractal character of 

Q ' the fluctuations and provides a new method to compute the localization length 

2 ' in terms of the fluctuations. In the case of a linear rotor, the fluctuations are 

CJ . independent of the kicking parameter k and exhibit self-similarity for certain 

values of the quasienergy. For given k, the asymptotic localization length is a 
I^ \ good characteristic of the localized lineshapes for all quasienergies. This is in 

Vh ' stark contrast to the quadratic rotor, where the fluctuations depend upon the 

strength of the kicking and exhibit local "resonances". These resonances re- 
sult in strong deviations of the localization length from the asymptotic value. 
The consequences are particularly pronounced when considering the time evo- 
lution of a packet made up of several quasienergy states. 
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Dynamical localization is an important manifestation of the quantum suppression of dif- 
fusive classical motion resulting from nonintegrable dynamicsETu. As the name suggests, the 
mechanism is analogous to the Anderson description of a low- dimensional, low-temperature 
insulator phase in terms of tight-binding models (TBM)I§Q. The relationship between these 
two seemingly disparate systems was made explicit^ in a class of kicked quantum Hamilto- 
nians of the form (in dimensionless units): 

H = K{p) + V{e)J26{t-n), (1) 

where K(p) denotes a general kinetic energy operator. Note that time is measured in units 
of spacing between kicks. The time-periodic nature of the Hamiltonian allows the time- 
dependent solution to be expressed, in terms of the one-step evolution operator U, as 

Uij{9,t)=ij{9,t + 1). (2) 

For kicked systems, U takes on the particularly simple form: 

U = exp {-iK{p)/n) exp {-iV{e)/h) , (3) 

where h refers to the effective quantization scale in dimensionless units. Further, the evo- 
lution of any initial condition can be expressed in terms of quasienergy states 0^ which 
satisfy 

U<P^ = e-'^ct>^ , (4) 

where the quasienergy uj is real as U is an unitary operator. 

The relationship between these quantum kicked systems and TBM becomes clear on 
projecting the quasienergy states onto eigenstates of K{p). In our context, these are angular 
momentum states. The equation satisfied by the projection coefficients Um onto the im}^ 
angular momentum state iso 

TmUm + Yl ^n-rUr = , (5) 

r 

where 

T^ = tan[(c^-ir(m))/2], (6) 

and the Wm are the Fourier weights of 

w{e) = -i^Ti{y{e)/2h) , (7) 

with respect to the angular momentum basis. This transformation provides a simple method 
to understand dynamical localization and recurrences in energy in kicked rotors. 

In this mapping, the integer angular momentum quantum number of the rotor corre- 
sponds to the lattice site in TBM. The free phase evolution between kicks provides the 
pseudorandom diagonal (on-site) potential while the kicking potential V{9) determines the 
range and strength of the hopping. Thus, under certain conditions, boundedness and re- 
currences in energy in the kicked rotor manifest themselves in quasienergy states which are 



exponentially localized on a lattice labelled by the angular momentum quantum number of 
the rotor. 

/^From a practical standpoint, this method of studying kicked rotors is particularly 
useful when the TBM contains only nearest-neighbor couplings. This has motivated 
many studie^B of a special class of kicked rotors where the potential is chosen to be 
V{6) = —2haictaii{kcos{6)), resulting in a TBM, 

k 

TmUm + 2 ("m+1 + ^m-l) = . (8) 

Note that we, unlike earlier treatments, explicitly retain the presence of the quantization 
scale h in the definition of the potential. This makes it clear why the classical limits of 
quantum rotors corresponding to this choice of potential are trivial. The TBM, however, 
are perfectly well-defined and serve as useful illustrative examples. 

This nearest-neighbor TBM has been studied^ for linear and quadratic rotors which are 
described by K{p) = ap, where a is an irrational number, and K{p) = p'^/2. In the diagonal 
term in the TBM, these translate into K{m) = am and K{m) = Tirn^ /2 respectively. Thus, 
h does not explicitly appear in the TBM analysis of the linear rotor. In the quadratic rotor 
we set, U = Sttct in keeping with the requirements of dynamical localization^. 

The linear rotor, where the diagonal disorder is quasiperiodic, was solved exactly. In 
particular, the density of states was shown to be identical to the average density of states 
for the Lloyd model of disordem for which the T^ are independent random variables with 
the Cauchy distribution 

P(Tn.) = -y^ , (9) 

7rl + T^ 

Furthermore, the numerically computed localization length was found to be in good agree- 
ment with the analytic, ensemble averaged, inverse localization length 7 of the Lloyd modelQ, 



cosh(7) = Vk-^ + 1 . (10) 

The localization length 7"^ depends upon k and controls the transient and recurrence time 
scales in the system. For the quadratic rotor where the diagonal disorder is pseudorandom, 
statistical studies of the sequences T^ showed similarities with the random potential with 
Lorentzian distribution. It was argued that in spite of some correlations, the quadratic 
model also exhibits the localization characteristics of the Lloyd modePa. However, our 
work here shows that these residual correlations have a profound impact on the fluctuations 
in the lineshape. These, in turn, can lead to strong deviations of 7"^ from the ideal Lloyd 
model prediction. 

In this paper, we use a novel method to directly compute the fluctuations in the ex- 
ponential lineshape for linear and quadratic rotors. We compare and contrast the linear, 
integrable, system with equally-spaced quasienergies, ujj = j *a mod(l) ( independent of k ) 
and the quadratic, nonintegrable case, where the quasienergies depend upon k. We reiterate 
that the classical limit of the Lloyd model is not well-defined, and so integrability or not 
refers only to the distribution of the energy levels. In other words, the terminology merely 
distinguishes between the two cases we consider and does not imply a classical limit. 



Our motivation is to understand better how the two different K{p) in the Icicked rotor are 
manifested in the locahzation properties of the equivalent TBMs. We characterize the differ- 
ences in these two rotor systems by using a recent techniquea for studying the fluctuations 
in the respective localized quasienergy states. We begin by factoring out the exponential 
envelope. Thus the projections of the eigenstates of the rotor Um are written as 

u^ = e-^l'"!//^ , (11) 

where rj^ are the fluctuations in the localized states at the m*^ lattice site. Thus, the 
fluctuations 1]^ can be related to the fluctuations in the asymptotic localization length 7, 
(where Um = e~^ ) as 

1-7 = log\Vm\/m (12) 

As explained below, we use an exact decimation scheme to compute the scaling properties 
of the fluctuations 77^ thereby establishing the fact that they are fractal. Furthermore, the 
scale length 7"^ need not be presupposed but can be self-consistently determined. It should 
be noted that this equation is valid for any TBM including random systems and therefore the 
method described above can be used to compute localization length for any nearest-neighbor 
TBM irrespective of the nature of diagonal or off-diagonal disorder. 
The fluctuations 1]^ satisfy the following TBM 

e~'^T]m+i + e'^r]m-i + XTmilm = , (13) 

where A = 2/k. We apply a recently developed decimation methods to this TBM with a in 
Tm was taken to be the inverse golden mean (v^ — l)/2. In this approach, the incommen- 
surability of the lattice was exploited by decimating all sites except those labelled by the 
Fibonacci numbers F„. This renormalization group approach was shown to an extremely 
useful tool to demonstrate self-similarity and to obtain universal characteristics of quasiperi- 
odic systems. Here, we apply this formalism to a linear rotor, described by quasiperiodic 
TBM, as well as to the quadratic rotor which is not quasiperiodic. We demonstrate that 
the Fibonacci decimation scheme is a very efficient method to compute fluctuations in the 
localization lengths irrespective of the nature of the aperiodicity of the TBM. Note that for 
non- quasiperiodic problems, Fibonacci decimation can be replaced by a more conventional 
one where every other site is decimated. 

After the n*'' decimation level, the nearest-neighbor TBM connecting the fluctuations at 
two neighboring Fibonacci sites is 

fn{m)r]{m + F„+i) = r]{m + F„) + en{m)r]{m). (14) 

The additive property of the Fibonacci numbers provide exact recursion relationdaor the 
decimation functions e„ and /„: 

f f^ _ fn-lii + Fn)fn{i + Fn) , . 

A = e„_i(z + F„) + fn-iii + Rn)en{i + F„). 



These can be iterated to machine precision as they do not depend on any parameter which 
could hmit the precision. It turns out that for the locahzed phase where 7 is always greater 
than zero, the decimation function /„ vanishes asymptotically and hence the resulting renor- 
malization flow simplifies to 



in+iyi') = -Cn-l 



[i + Fn)en{i) . (17) 



In view of this further simplification, the above equation can be iterated up to 35 iterations 
which corresponds to studying TBM of size up to 14,930,352. 

Any fractal character in the fluctuations can be inferred by non-trivial asymptotic behav- 
ior of the functions e„. In particular, the convergence of the renormalization flow either to 
a non-trivial limit cycle ( which implies self-similarity ) or a strange set is a clear indication 
of fractality. 

It is interesting to note that the decimation function e„ also determines the localization 
length as 7 — 7 = log |e„|/F„. When the 7 is not known, the above equation determines it 
self-consistently to a very high precision. 

In order to calculate the exponential lineshape for a given quasienergy , we iterate the 
TBM (Eqn.(8)). Starting from the site at +A^, we iterate the equation backwards and 
simultaneously iterate forwards from site — A^. We then match these backward and forward 
iterates at m = by adjusting the phase factor u, thereby determining the quasienergy. 
Note that in this method, the localization center is always at ?7i = 0. 

An extremely accurate method to determine the quasienergies results from rewriting the 
TBM as a quasiperiodically driven mapliil. This is obtained by defining 

Xm = Um-l/Um , (18) 

which transforms the TBM to 

_ -1 

''™+' ~ x^ + Atan((cu-K(m))/2) ' ^^^' 

The localized phase of the TBM manifests itself as a strange nonchaotic attactor of this map, 
reflected in the expression for the Lyapunov exponent of the map jj, = — 27^11. The negative 
sign is crucial as it implies that an attractor with /x < corresponds to diverging Um with 
increasing m. Similarly, an attractor of the inverse map corresponds to diverging Um with 
decreasing m. Thus, an exponentially localized function is obtained by starting from the 
intersection points of these two attractors and reversing the direction of iteration in each 
case. The quasienergies of the kicked rotor are used as tunable phase factors to ensure that 
the intersection points of the attractor always occur at ?n = 0. This provides an extremely 
accurate method to compute the quasienergies. The resulting u^ is exponentially localized 
in both directions with the localization center ttiq always at zero. Note that a consequence 
for the results we show is that the lattice site label is always relative to the localization 
center. 

As stated earlier, the primary difference between the linear and quadratic rotors lies in the 
character of the fluctuations. In both cases, the rotor wave function Um, with exponentially 
decaying envelope, exhibits fractal fluctuations r/^ which decay as a power law, rj^ ~ rn^ . 
The exponent p is related to the decimation function as p{n) = ln[abs{en)]/ln{Fn). Asymp- 
totically ( in the limit of large decimation level ), e{n) ^ v{Fn), the exponent p is a measure 



of the fluctuations in the exponent 7 as p{n) = (7^ — ^){Fn)/ln{Fn). Note that the mea- 
sure p{n) depends on the decimation level n. The behavior of p with kicking parameter k 
distinguishes the two rotor cases we consider. 

Figure |l| shows the self-similar fluctuations 77^, for tu = 0, in the case of the linear rotor. 
It should be noted that the fluctuations are described by a non-trivial period-6 limit cycle 
of the renormalization flow, that is rjp,^ = i]F„+e ■ Further, the same period-6 was found for 
all values of k thereby establishing the fact that the fractal fluctuations are independent of 
k for k > 0. The multifractal nature of these fluctuations was also confirmed by the /(a) 
curve. For other values of quasienergies, the fluctuations are not described by a limit cycle. 
Instead, the decimation functions e„ (for all quasienergies) converge on an invariant set with 
fractal measure, independent of k. For our purposes, it is more illustrative to focus on the 
impact of these features on the localized lineshape. Given the relationships constructed 
earlier, this universality implies that the convergence of the localization length 7"^ to the 
asymptotic value 7"^ is also independent of k. 

As seen from the Fig. 0(a), the universality in the linear rotor is reflected in p{n) which is 
independent of k at different decimation levels. In contrast to the linear rotor, the quadratic 
rotor shows fluctuations which depend strongly on k. As seen from Fig. 0(b), the power p{n) 
exhibits many spikes as a function of the parameter k. These spikes are the result of local 
"resonances" in the fluctuations. The nature of these resonances is seen from Fig. ^(c) where 
the magnitude oir]^ changes by several orders of magnitudes within a few sites. The location 
in m and the magnitude of these 'jumps' depend on both k and the quasienergy u. Two 
values of k are shown in Fig 0(c) with the dark curve {k = 3.252) corresponding to a larger 
spike (for both n) in panel (b). Figure ^(d) shows that the resonances lead to 'shoulders' in 
the exponential lineshape, indicating local variation in the exponential envelope. It is worth 
noting that the quasienergies vary with k and are computed for each k using the matching 
condition described earlier. 

The location of the resonance determines the impact of these local deviations. If the 
resonance occurs close to the localization center, then the deviations in 7 from 7 are clearly 
significant. However, if the resonance site is far from the localization center, then the expo- 
nential envelope diminishes the importance of the deviations. Figure ^ shows the lineshapes 
associated with single quasienergy states for both rotors. As seen from the linear rotor 
lineshape in Fig. |5(c), the absence of resonances means that the Lloyd model estimate is 
a very good approximation to the actual calculated localization length, all the way to the 
localization center. This result was found to be true for all quasienergies and for all values 
of the kicking parameter k. However, for the quadratic rotor shown in Fig. ^(a) and (b), 
the Lloyd model estimate of 7"^ appears to be correct only in the tails of the localized line 
shape for most of the quasienergies. Deviations from the asymptotic 7"^ depend strongly 
on both the parameter k and the quasienergy, as seen by contrasting panel (a) and (b). For 
u = 0.17, 7~^ is clearly a better fit for the associated quasienergy state than for case of 
cu = 0. 

Therefore, localization in the quadratic rotor starting from an initial wave packet com- 
posed of several quasienergy states should not be well described by the Lloyd model. To 
illustrate this, we consider the evolution of a plane wave (at m = 0) under the repeated 
action of the single-step evolution operator U. The canonical method of Fast Fourier Trans- 
forms is used to get the localized lineshape after a large number of kicks. The probability 



distribution /{m) across lattice sites m is then constructed. As seen from Fig ^, the hnear 
rotor ((b) and (d)) coincides with the Lloyd model prediction all the way to the localization 
center. This is in stark contrast to the case of quadratic rotor in panels (a) and (c) where the 
deviations are strongest close to the localization center. It is also evident that the magnitude 
of these deviations depend on k. Exponential fits near the centers of the lineshapes yield 
values of 7"^ which are clearly different from the asymptotic estimates given by 7"^. 

The conditions under which the Lloyd model calculationEi was made allow us to speculate 
on the possible reasons for both the resonances and the associated deviation in 7. The 
calculation of the estimate 7 required the T„ to satisfy a specific distribution. However, 
dynamical phase correlations would lead to a violation of this requirement. This was verified 
both by constructing the return mapping for the on-site potentials and by directly plotting a 
histogram of on-site terms and contrasting with the required distribution. These were noted 
in earlier work as wellBQ and are clearly different in the linear and quadratic cases. 

We believe that the local resonances and fluctuations seen in the localization character- 
istics of the quadratic rotor are generic to pseurorandom systems. We have verified that 
TBMs with bounded onsite potentials such as cos{2TTamy) also exhibit the characteristics 
similar to that of the quadratic rotor for v > 2. u = 1 constitutes a special case as the model 
reduces to the well-known quasiperiodic Harper equationEJ where the fractal localization 
character was also found to be independent of the couplingp. It should be noted that the 
Harper equation has been recently solved using the Bethe-AnsatzEd implying some sort of 
'integrability' in the model. In view of this, we speculate that for the linear rotor, the k 
independence results from the integrable nature of the problem. 

Ideally, in quadratic rotors, the possibility of correlations can be recognized from study- 
ing the classical dynamics resulting from the kicked rotor Hamiltonian. Specifically, in the 
context of quantum dynamics in mixed phase spaces, quantum phase correlations can be 
associated with the presence of invariant structures in the classical phase space. This rela- 
tionship is of great current interest in the new context of quantum manifestations of classical 
anomalous transportlljO. Recent work suggests that there may be a relationship between 
the large fluctuations in the localization length and the anomalous diffusion in the classical 
phase spaceo. The lineshapes also exhibit shoulders similar to the ones shown here. We pro- 
pose to examine more closely the association of quantum phase correlations with structures 
in the associated classical phase space. However, as mentioned earlier, this is not possible 
in the special class of rotor studied here as the corresponding classical mapping is not well 
defined. In keeping with this general motivation, we are presently extending our work to 
models where this is not an issue. 
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FIGURES 

FIG. 1. Self-similar fluctuations in a linear rotor for uj = . The fluctuations repeat themselves 
at every G*'' Fibonacci site, that is -F(6),-F(12), • • • etc.. This type of "translational invariance" in 
Fibonacci space is described by a period-6 limit cycle of the renormalization flow. Note that this 
behavior is independent of the kicking parameter k. 

FIG. 2. The variation in the power law exponent p{n) as a function of the kicking parameter 
k is shown for (a) the linear and (b) quadratic rotors at two different decimation levels. The darker 
and lighter curves respectively correspond to 12 and 15 decimation level. Panel (c) shows the 
fluctuations rjm versus m for the quadratic rotor at A; = 3.252 (darker curve) and at k = 2.9 (lighter 
curve). Note that the larger value of k corresponds to a peak (spike) in (b) while the other does 
not. The jump in magnitude of the fluctuations is seen in the corresponding lineshapes shown in 
(d), with larger fluctuations leading to 'shoulders' seen in the lineshape. These 'shoulders' result 
in variations in localization length from the corresponding asymptotic values. 

FIG. 3. Single localized quasienergy states (points) for the quadratic rotor (a) and (b) and the 
linear rotor (c). The parameter k = 2.8. Note that (a) and (b) correspond to quasienergies 0.00 
and 0.17 respectively. The linear rotor case is computed for to = 0. The Lloyd model prediction 
for the localization length 7"^ is plotted in each case (solid line) to guide the eye. 

FIG. 4. Localized lineshapes (points) starting from a plane-wave initial condition for linear 
(panels (b) and (d)) and quadratic rotors ((a) and (c)). The first column considered k = 2.8 while 
the second is for k = 5. In all cases, the solid line indicates the Lloyd model result. 
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